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Abstract 



A theoretical study of vesicles of topological genus zero is presented. The 
■ bilayer membranes forming the vesicles have various degrees of intrinsic 

ON ■ (tangent-plane) orientational order, ranging from smectic to hexatic, frus- 



trated by curvature and topology. The field-theoretical model for these 're- 



cti 

e 

*^ • atic' surfaces has been studied before in the low temperature (mean-field) 

c . 

limit. Work presented here includes the effects of thermal fluctuations. Us- 



ing the lowest Landau level approximation, the coupling between order and 
shape is cast in a simple form, facilitating insights into the behaviour of 
vesicles. The order parameter contains vortices, whose effective interaction 
potential is found, and renormalized by membrane fluctuations. The shape 
of the phase space has a counter-intuitive influence on this potential. A cri- 
terion is established whereby a vesicle of finite rigidity may be burst by its 
own in-plane order, and an analogy is drawn with flux exclusion from a type-I 
sup er conductor . 

PACS N%: 82.70, 02.40, 68.15, 64.70.M 
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I. INTRODUCTION 



When amphiphilic molecules (those with hydrophobic and hydrophilic parts) are dis- 
solved in water, they hide their aliphatic tails by grouping together and thus assemble 
themselves into structures. Depending on the geometry and chemistry of the particular 
surfactant, these structures may be anything from nanometre-scale mycils to macroscopic 
membranes arranged in stacks, bicontinuous networks of pipes, or closed surfaces. The 
literature on spontaneously self-assembling amphiphiles and biological and liquid- crystal 
membranes is vast. See for example |JT|, M or |3|]. 

Research presented here is concerned with bilayer fluid membranes with internal ori- 
entational degrees of freedom. Fluidity of the membrane refers to the fact that molecules 
can move freely within the surface, although perpendicular deformations have an energy 
cost. The molecules have no positional crystalline order as the temperature is above the 
crystalline-to-fluid transition at which lattice dislocations proliferate. However, for some 
liquid crystals, this is not coincident with loss of orientational order, which occurs at a 
higher temperature. Various kinds of orientational order are possible. Liquid crystals in the 
smectic-C phase have carbon-chain tails which are tilted with respect to the local normal to 
the surface. The local mean direction of tilt, projected onto the local tangent plane of the 
membrane, gives a two-component vector order parameter, which disappears at the continu- 
ous transition to smectic-A phase, where rotational symmetry is restored . Alternatively, 
the orientational order may be hexatic. In this case, the Kosterlitz-Thouless (K-T) tran- 
sition at which lattice disclinations proliferate, destroying the quasi-long-range hexagonal 
bond orientational order, is non-coincident with the lower temperature K-T transition at 
which dislocations proliferate ||. Nearest neighbour bond orientations give a basis for defin- 
ing an order parameter, which again has two components in the local tangent plane. But a 
six-fold ambiguity in the choice of nearest neighbours arises from the local six-fold rotational 
invariance of the membrane. Mapping the local tangent plane onto the Argand plane, with 
some arbitrary direction chosen for the real axis, allows the two-component order parameter 
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to be well defined clS £L complex number thus: 



V(o-) = (exp(i6e(ff))) , 



where cr is some two-component coordinate defined on the surface, and G is the bond angle 
with respect to an arbitrary coordinate axis. Notice that this definition is independent of 
which of the six bonds is chosen, since rotations of the membrane through 7r/3 correspond to 
a phase change of 2tt in ip, or an identity transformation. Hence, this two-component order 
parameter has the appropriate six-fold local rotational invariance. A similar complex order 
parameter can be defined for the smectic-C order described above. This has local full-turn 
invariance only, and magnitude tJ)q. Hence it may be mapped onto a complex field thus: 



where G gives the orientation of the tail vectors projected onto the tangent plane. In each 
case, () denotes a local thermal average. The notion may be generalized to an in-plane 
orientational order parameter with local ra-fold rotational invariance being represented by a 
complex field 



where G gives the orientation of one of the n principal directions of order. For smectic 
and bond-angle order in liquid crystals, n G {1,2,4,6}. Geometrical arguments for this 
restriction on the possible orders of rotational symmetry are given in |J. 

On a flat membrane, the continuous development of ip is well modelled by a Ginzburg- 
Landau free energy functional of the form 



where dA is an element of area, and r, u and C are parameters of the particular conditions 
of the system being modelled. 

On a flexible membrane however, the free energy functional is a little more complicated. 
Let the membrane be treated as an ideal, two-dimensional, smooth surface. A Hamiltonian 



-0c (tr) = (^ exp(iG(cr))) , 



if)[tr) = (^ exp(mG(cr))) , 
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due to Helfrich |7j governs elastic deformations of an isotropic (ip = 0) fluid membrane. 
This Hamiltonian is a functional of the intrinsic and extrinsic curvatures of the membrane, 
as defined below. Let R(cr) be the position vector in R 3 of a point on the surface with 
two-dimensional coordinate cr = (o"i,o"2). The surface has a metric tensor with components 
9ab — d a R ■ df,R where a and b label the coordinates. The metric has inverse g ab and 
determinant g. A covariant curvature tensor is defined by K a b = N ■ d a dbR where N is the 
local unit normal. Note that dot products are evaluated in R 3 . Using Einstein summation 
convention, and the metric and its inverse to lower and raise indices in the standard way, 
a more useful curvature tensor is defined thus: K% = g ac K c i,. Its determinant K is the 
Gaussian curvature of the surface, also called the intrinsic curvature, as it is experienced 
by 'flatlanders' and other physical entities living within the two-dimensional world of the 
surface, who are oblivious to the way in which the surface is embedded in 3-space. Its 
trace K% however is an extrinsic property, which cannot be determined by flat lander s[]. This 
quantity is known as the total curvature, or twice the mean curvature, and is the sum of the 
two principal radii of curvature of the surface at a given point || . Helfrich's Hamiltonian, 
describing the elastic properties of a constant area fluid membrane, can now be expressed 
as 

^Helfrich = / {HKf + K G K}^g-d 2 a. (1.3) 

Note that ^Jgd 2 a = dA. The phenomenological constants in this expression are the bending 
rigidity k and the Gaussian curvature modulus kq. On a closed surface, the second term in 
equation (|1.3j ) is a topological invariant, according to the Gauss-Bonnet formula 

J KdA = 2n X 

where \ is the Euler number of the surface |10| . For an orientable surface of genus (number 
of 'handles') G, the Euler number is given by 

x The idea of 'flatlanders' was used by E. A. Abbott in his book "Flatland" ||. They are hypo- 
thetical beings who live in a two-dimensional space, and have no conception of a third dimension. 
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X = 2(l-G). 

Hence the dynamics of a membrane of given topology are not influenced by the second term 
of equation ( |1.3|) . What follows will concentrate on membranes of spherical topology (genus 
zero), so this term will be dropped. 

The free energy functional governing the dynamics of n-atic fluid membranes cannot be 
found by simply adding together equations (|1.2|) and ( |1.3|) , since ( |1.2|) contains a gradient 



operator, which must be expressed in terms of coordinates which are no longer flat, but 
curvilinear. The derivatives in the flat-membrane expression | V^| 2 = d x ip*d x ifj + d y ifj*dyip 
must be replaced by covariant derivatives thus: |V"?/>| 2 = D* a ip*D a il) where D a is of the form 
d a — inA a with the connection A defined such that the free energy is invariant under gen- 
eral coordinate transformations. Clearly, A itself is not coordinate-invariant, since ip has 
coordinate dependence due to the arbitrariness of the reference axis for B in equation ( p..l|) . 
There is some freedom in the choice of the gauge field A. Here, A is chosen to be the 'spin 
connection' A a = ei ■ d a e 2 where is a unit vector in the direction of d^R- Finally, we have 
the free energy functional for a closed, n-atic fluid membrane: 

F[i/>{<r),R(tr)] = J d 2 a^{r\^\ 2 + \u\^ + CD^D^ + \^K a a ) 2 } (1.4) 

and, since the membrane's total area is fixed at A by the number of incompressibly packed 
constituent molecules, there is also a constraint 

J d 2 a^ = A. (1.5) 

The covariant derivatives provide a coupling between order and the shape of the membrane. 

So far, the model and formalism are identical to that used by Park, Lubensky and 
MacKintosh in their mean-field treatment of n-atic fluid membranes of genus- zero (vesicles) 
||1 1|| . They noted that, not only is n-atic ordering partly frustrated by curvature of its two- 
dimensional space, making parallelism impossible to achieve, but it is further frustrated by 
the spherical boundary conditions. Orientational order intrinsic to a surface of Euler number 
X must have topological defects whose indices sum to x- F° r instance, smectic-C (vector) 



order can form defects (sometimes called vortices) of index 1. Hence, on a sphere (Euler 
number 2) and in the absence of any anti- vortices (index -1 defects), there must be two 
defects in a Sm-C order parameter. At these poles, tp vanishes smoothly to avoid infinite 
gradients. Clearly, an n-atic can form defects of index 1/n since it can rotate by less than 
a full turn in circling a pole, without its phase slipping. So n-atics form 2n vortices on a 
sphere, in the absence of excited vortex-anti-vortex pairs. 

A similar approach to that of Park et al. will be used in this paper. The small defor- 
mations from sphericity of the n-atic fluid membrane are expressed as a real, scalar field, 
and the order parameter as a complex field ijj as described above. The shape deformation 
field is expanded in spherical harmonics, and ip is expanded in eigenf unctions of the gradient 
operator. As an approximation, the series expansion for ip is truncated and only the de- 
generate set of lowest-eigenvalue functions are kept. The same approximation was used by 



Landau in finding the wavefunction of an electron in a magnetic field |12[], and by Abrikosov 
in his treatment of flux lattices in superconductors [P~5 | . The study diverges from that of 
Park et al. when a simple expression is calculated for the coupling between the two fields. 
The free energy is cast into a form which allows the shape fluctuations to be integrated out 
exactly, giving renormalized coefficients for the order parameter. This effectively maps the 
deformable sphere problem onto a rigid sphere problem, which is soluble for the n = 1 case. 
In section [ED], the stability of the system is analyzed to ascertain the conditions under which 
topological defects become too expensive, and are excluded from the membrane in a similar 
fashion to flux exclusion from a type-I superconductor. The partition function and various 
expectation values for n = 1 are calculated in section [V]], using some special symmetries 
of the problem. For other values of n, the system is not soluble, but approximate expecta- 
tion values are found in section |V11| , using diagrammatic expansion and the Hartree-Fock 
method, which gives correct results in the limit of high temperature. 
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II. CALCULATION 



Using the expression D a = d a — inA a , it is easy to show that, on a closed surface, 

Both sides are real for all ip hence the operator —g^^D a {^fgD a ) is Hermitian. Eigenvalues of 
this operator scale as A. A useful set of orthonormal basis functions is formed by normalized 
eigenfunctions of this operator on the unit sphere, where y/gd 2 a = dfl which is the usual solid 
angle element (sinOdOdip in spherical polar coordinates). These functions are sometimes 
referred to as 'Landau levels'. There are (2n + 1) degenerate Landau levels with the lowest 
eigenvalue, n. If ijj is expanded in this complete set of functions then, close to the mean-field 
phase transition, the partition function is dominated by configurations involving only the 
lowest Landau levels. Hence, in this regime it is a good approximation to expand ip in 
these complex, degenerate lowest levels only. Thus the number of degrees of freedom in ijj 
is reduced to 2(2n + 1), corresponding to freedom in the positions of the 2n vortices in the 
two-dimensional surface, together with an overall complex amplitude. On the unit sphere, 
D ( /aD a ) 

avyy — )_ = _f d0 Q e + cosec 2# q q + cot e q $ + 2 i n cosec q co t B d v - v? cot 2 9). (2.1) 

The lowest eigenfunctions of this operator and similar operators for other gauges, are pre- 



sented in different forms in Jl4]], [ffH] and [11], and are re-expressed here as 



U " ' 1)1 sin n+p (^)cos n - p Aexp( W ) (2.2) 



Yp \ An(n + p)\(n-p)\ y 2 J y 2' 

for integer values of p between — n and n. Any function ip which is a linear combination of 
these functions, obeys the relation 

dgip = (ncot9 — i cosec 9 8^) ip (2.3) 

Let deviations from sphericity of the vesicle's shape be parameterized by a real, dimen- 
sionless, scalar field p{cr) in the following way. If the ground-state vesicle shape is a sphere 
of radius Rq then, with the origin at the centre of the vesicle, 



\R(*) \ = (1 + p(<r)) R . (2.4) 
This is a 'normal gauge' parameterization, and therefore carries the correct weight in a 



statistical ensemble [16 



It will emerge that < p 2 >~< |^| 4 >. Equation ( 1.4|) approximates the potential by a 
series expansion, truncated at fourth order. Hence, it is consistent to also truncate to order 
p 2 and \ip\ 2 p. ie. the vesicles under consideration deviate little from spheres. To this order, 
y/g\i/)\ 2 = R 2 sm9(l + 2p)\^\ 2 and 

—tp*D a ( y /gD a )'ijj = n\ip\ 2 sm9 + 2nip* {id v p dgip — idepd v ijj + 2ncos9vp dgp}. 

Applying the relation Q2.3) and integrating by parts gives 

J -rD a (^D a )^d 2 a = nj d^\ 2 {\ - V\p) + 0{\^\ 2 p 2 ). 



where V 2 = (dgdg + cosec 2 # d^d^ + cot 9 dg) is the covariant Laplacian. So the expression 
for the free energy becomes 

F = R 2 Jd{l{(r- r c )H 2 (l + 2p) + \u\^\ 4 + r,#| 2 (2 + V 2 ± )p} + tt He lfrich + & 
where 



c — p2 ' 



Cn 

W 
if) 

It also emerges that \ip\ 2 ~ (r — r c ) below the mean-field transition. See section for more 
explanation. Furthermore, from the definitions of K% and p, it follows that 

K = - -J-(i - P - \v\ P + P 2 + pV\ P ) 



r 







and, from equation (|1.5|) . 



A = A<kRI + R 2 J dQp(l + \V\)p 



since p and Rq are defined in equation (|2.4[ ) in such a way that J pdQ = 0. 
Hence, to C( I "0 1 4 ) ? we ma Y write 



AimC 
rc = ~ — 

and the free energy 

Fty, p] = j dn {(r - r c )|^| 2 + i M |^| 4 } - Cr#| 2 (2 + Vi)p + « {2 + pVip + |(Vip) 2 } 

(2.5) 

with the field p unaffected by the fixed area constraint to this order. 

The coupling between the fields ip and p, describing order and shape, has been cast into 
a very simple form, using some special properties of the lowest Landau levels of tp in the 
spherical geometry. Note that a similar restriction has not been put on the phase space of p. 
All modes of shape fluctuation are available for small amplitude excitation. The remarkable 
simplicity of equation ( |2.5|) is conducive to further exploration of the properties of n-atic 
vesicles, hitherto hindered by unwieldy notations. 



It is clear from equation ( |2.5| ) that the mean-field transition temperature has been shifted 
from r = to the lower temperature r = r c due to the non-zero lowest eigenvalue of the 
operator in equation ( |2.1|) in the spherical geometry. Hence, ordering is frustrated by the 
curvature of the space. 

Now let p be expanded in eigenf unctions of the covariant Laplacian Vj_. These functions 
are spherical harmonics Y™, defined by 

V\Yl n = -1(1 + 

The s-wave (I = 0) harmonic is not included in the series expansion, since this mode of 
deformation is already represented by re-scaling Rq in equation (|2.4j ). The I = 1 modes are 
also excluded, as F will turn out to be independent of them. These three modes describe a 
positive deformation on one hemisphere of the vesicle, and negative on the opposite hemi- 
sphere. Hence, they simply correspond to the three degrees of translational freedom. The 

complex coefficients pi m become the dynamical variables, where 

00 1 

1=2 m=-l 



The deformation field p is constrained to be real by demanding that 

PL = {-^Pl-m 

which follows from the identity Y{ a * = (-l) m Yf m . 

As stated earlier, the order parameter field ip is expanded in the basis functions given in 
equation (|2.2|) thus: 

[4tt n 

^ = \ ~A H a P^P ( 2 ' 6 ) 
» 1/1 p =-n 

where a p are (2n + 1) dimensionless, complex, dynamical variables. 

The free energy can now be expressed in terms of the dynamical variables thus: 

F = a a* a p + ^ P pqrs a*a*a r a s + ^ l m imP^fl q Pim + Y ^iP*i m Pim (2.7) 

P pqrs pqlm Im 

where 

a = r + — — (2.8a) 



P pqrs = —r \ <f>;<f> q <f>r<f>s dn (2.8b) 



2nu 
~A 

l pq im = + m -I) J ^ q Y™ dn (2.8c) 

A/ = \kI{1 2 - + 2) (2.8d) 

all of which coefficients are real. Indices run over the intervals — n < p,q,r,s < n and 
2 < / < 00 and —l<m<l. Notice that, as stated earlier, I = 1 spherical harmonics 
contribute to neither the curvature energy nor the coupling, due to the factors of (I — 1). 
Notice also that equation ( |2.7[ ) is quadratic in the variables pi m . Hence, in a partition 
function defined by 



Z = J e~ F « a ^^» [] da* q da q dp* Vm ,d PVn , 

q I'm 1 



the integral over all shapes may be done explicitly by completing the square, using the 
identity •y pq i m = (—^) m lq P i-m- Hence a new effective free energy, defined by 
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Z = f e- Feff « a *»l[da* q da q (2.9) 



is found to be 



F eff = «E a> p + £ ^ rs a>;a r a s + ft p (2.10) 



p pqrs 

where 

1 OO £ 

and Hp is the free energy of shape fluctuations due to the Helfrich Hamiltonian alone, given 
by 



1=2 



1 / A « 



2 4 n 



Ci In k + C2 



where (1 and (2 are infinite constants in the continuum limit. This infinite part of the free 
energy will henceforth be renormalized away. 

Consider for a moment the coefficients 7 p(? /m, defined via the integral / (p^gY^dQ. Now, 
the spherical harmonic Y^iO, <p) is of the form P[ n (cos9) exp(irmp), where P™ is an associ- 
ated Legendre function, of the form sin m 9 J2k=Hi bk cos k 9 with constant coefficients The 
associated Legendre functions form a complete, orthogonal set of functions for each value of 
m. Note also that (fi*4> q oc sm p ~ q 9(1 - cos9) n+g (l + cos 9) n ~ p exp i(q — p)(p. It follows that 
it is possible to write the function as a finite sum of spherical harmonics, ie. with some 
constant coefficients c^, it can be written 4>*(p q = Y^jZo c jYf~ q * ■ Hence (j) p 4> g is orthogonal 
to all spherical harmonics not included in this sum. So •y pq im vanishes for I > 2n, and the 



infinite sum in equation ( 2.11 ) can be replaced by a finite sum 



/5&. = /W.-zE E (2.12) 

4 1=2 m=-l 1X1 

Physically, this means that only the modes of deformation between / = 2 and I = 2n couple 
to the lowest Landau levels of order. 
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III. STABILITY ANALYSIS 



The original system of orientational order on a fluctuating sphere has now been mapped 
onto a rigid sphere problem, in which a p are the only dynamical variables, by replacing the 
coupling j3 pqrs by the effective coupling j3 p g rs , renormalized by shape fluctuations. The free 
energy closely resembles that of a superconducting film, penetrated normally by magnetic 
flux quanta, around which the supercurrent flows in vortices. It is well known |l7j] that 



changing the sign of the fourth-order term from positive to negative changes the nature 
of the superconductor from type-II to type-I, from which all flux vortices are excluded. 
The corresponding phenomenon of excluding all vortices from an n-atic vesicle appears 
impossible, since the vortices are topologically imperative. There is, in fact, an analogous 
exclusion of vortices, brought about by a catastrophic change in the vesicle's topology, at 
which point the model breaks down, since the order and deformation fields become large. 
In the analysis of type-II superconductors, there is a similar break-down of the model in the 
marginal limit. The criterion for a 'type-I' vesicle is now found by a stability analysis. 



The free energy of equation ( |2.10|) is a function of the (2n + 1) complex variables a p . 
It describes a stable system if F is large at infinity in all directions in the phase space, ie. 
if F — > +00 as any linear combination of the variables a p goes to infinity. Such a path 
to infinity may be parameterized by a p = rjpE, where H — > 00. The constants r] p define a 
phase-space direction and may be kept finite, without loss of generality, by the constraint 



J2 P V P Vp = 1- Substituting into equation fl2.10| ) and neglecting quantities of less than fourth 



order in H yields the criterion for stability 



pqrs 



There is one particular set of constants r] p = 77° which minimizes this sum. This set definitely 
exists, due to the introduction of the constraint above. So the stability criterion may be 
expressed in terms of this specific set of constants thus: 



pqrs 
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But this set of constants which minimizes the sum, subject to an overall magnitude con- 
straint, is simply the set of mean-field values of a p , multiplied by some scale factor. Hence 
the criterion for stability is J2 f3pq rs ( )m.f. > 0, and may be written as 

/ iv^.dn > ^E (< t(?+i) 1} 1/ l^l-.^f- 

Let 'malleability' p be defined as the quantity n 2 C 2 /AuK, which is a measure of the ease 
with which the order may deform the vesicle. The marginal case, defining the crossover from 
'type-II' to 'type-I' behaviour is given by p = p c where 

2 



\l,xrdn 



Using symmetry arguments to establish the mean-field configuration, p c (n) is found for the 
first two cases to be: /i c (l) = 9/47T, // c (2) = 125/15677T. For vesicles this malleable, the 
order parameter can deform the shape from a sphere to an intrinsically flat cylinder or 
polyhedron, and hence ordering is no longer frustrated by curvature. In other words, the 
vesicle bursts. This interpretation may not be the whole story, since the model breaks down 
as this instability is approached, due to the large magnitudes of the fields ip an d p- But it is 
clear that 'ordinary' or 'type-II' behaviour is exhibited when p < p c and that a qualitatively 
different behaviour, requiring a different model, exists in vesicles whose malleability is greater 
than this finite threshold. It will be demonstrated later that ordering increases and the region 
of temperature over which it arises sharpens with increasing malleability. The relationship 
between the magnitude of this order and the degree of shape deformation is now calculated. 

IV. SHAPE CORRELATIONS 



The original model for genus-zero vesicles with n-atic order, used by Park et. al. ||TT 
embodied in equations (|1.4J) and ( |1.5| ), has been reduced to the much simpler form of equation 
fl2.5|) and all of the shape deformation degrees of freedom have been integrated out, leaving 



the effective free energy in equation ( 2.10|) , with coefficients calculable from equations (2.12) 
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and ( p.8| ). Park et. al. solved the model in the zero temperature limit (without fluctuations), 
and thus were able to produce diagrams of the vesicle shapes, from the fixed values of pi m 
for each value of n. In reality of course, no such fixed values exist, as the shapes and vortex 
positions are in a constant state of flux, but expectation values of various moments of the 
variables pi m and a p can be calculated. The variables pi m have been eliminated from the 
problem, so transformation equations must be found, relating their expectation values to 
those of a p . This is done with the use of equation ( |2.7|) , and integrals of the form 



< Pim >= \\ V[i!,p}pi m e~ 



with the results 



< Plm >= Ipqlm < 0* p d q > (4.1) 

pq 



and 



< Pl 1 m 1 Phm 2 >= -j T A 2^ Iqrhrmlsphm* < O p d d r d s > H (4.2) 

4A h l\ h pqrs I\ h 



where 5 xy is the Kronecker delta. The expression for the second moment (equation (|4.2|) ) 
is in two parts. The first represents deformations of the vesicle due to interaction with the 
order (ie. ip is attempting to align, and expel Gaussian curvature, except at the vortices 
where it is small). This part vanishes for higher spherical harmonics (I > 2n), which do not 
couple to the lowest landau levels. The second part represents uncorrelated deformations 
due to thermal excitation. 



V. VALIDITY OF THE APPROXIMATIONS 

Before proceeding further, an assessment is made of the regimes of validity of the ap- 
proximations employed so far. 

The Ginzburg-Landau model carries with it an implicit approximation. It contains a 
truncated series expansion of the potential acting on ip. For this to be valid, it is required 
that higher order terms in the expansion are smaller than lower order terms. This condition 
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is satisfied at high temperature where ip is small. It is also true at low temperature, provided 
that (/i c — ^i)>Tt~ 1 . Clearly the model is good far from 'critical malleability', as stated above. 

In using the lowest-Landau-level approximation, it is assumed that the 'ground-state' 
modes dominate the order parameter field. Let the expectation value of their coefficients 
be denoted (a* o a po ^ and that of the next lowest level be denoted (cC pi a P ^j- Then the ratio 
(a* pi a p -^j I (a* o a Po ^ must be much less than unity. An order-of- magnitude estimate of this 
quantity is found from a high-temperature calculation to be ~ a /(a — r c ), from which it 
follows that | a | <C C / A. So this approximation is valid in the neighbourhood of the mean- 
field transition. Let its region of validity cover a large domain in figures 1-4, so that all work 
presented here is correct. This requires that C 2 /Au ^> (/i c — //). 

One final approximation must be justified. In deriving equation Q2.5Q , a factor (1 + 2p) 
was dropped from the line above, with the explanation that \ip\ 2 ~ (r — r c ), so this term was 
of too high an order in small quantities. This is only apparent at low temperature, but the 
derivation holds true at all temperatures for the following reasons. Replacing this lost factor 
is equivalent to multiplying TpgZm 

by l + 2(l-^)/(7 + 2)(/-l). Demanding that this is close 
to unity leads to the condition \a\ <C C/A, which is identical to the lowest-Landau-level 
criterion. 

It remains only to calculate expectation values of the order-parameter field using equa- 
tion (|2.10| ). For n = 1 (vector order) this is done in section [Vj by making use of various 



symmetries of the system. For other values of n, the calculation cannot be performed ex- 



actly, so a further approximation will be introduced in section VII, where the problem is 



solved for general n. Either of these sections may be read without reference to the other. 

VI. NON-PERTURBATIVE SOLUTION FOR VECTOR ORDER 

The n = 1 case (Sm-C order), solved at mean-field level in ||, has just two vortices on 
the spherical surface, and therefore has more symmetries than higher-??, cases, which can be 
exploited in the solution of the integrals of equation (|2""I)|). 
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A vesicle with vector order has two topologically required defects in the order parameter 
field. It is modelled using three complex dynamical variables, giving six degrees of freedom, 
corresponding to the positions of the two vortices in the two-dimensional space, plus an 
overall complex amplitude. Its energy is invariant under four obvious symmetries. The three 
Euler angles which relate to the spatial orientation of the vesicle are clearly irrelevant degrees 
of freedom, as is the complex phase of the overall amplitude of the order parameter field, 
which describes a global rotation of the orientation of in-plane order. The two remaining 
relevant dynamical quantities are the real amplitude of the order parameter and the geodesic 
distance between the vortices. A transformation of variables is now performed on equations 



^9|) and ( 2.10|) to make these symmetries explicit. 



Let us re-express the lowest-Landau-level order in the form used in |TT 

$ = , /lL II f sin - cos 6 ±e^-^' 2 - cos - sin e ± e -<v-^)A 
V A k=1 \ 2 2 22 J 

where (6k, (pk) is the position of the k th vortex in spherical polar coordinates. The factor 
^/(12/A) is convenient for the definition of the overall complex amplitude if> . Equating this 
expression to equation (|2.6|) produces the transformation equations: 



Ql = 2^ COS - COS h e ~^+^)n ( 6-la ) 

a_! = 2^ sin j sin d -l e ^+^)/2 
a = ^ 



V2^ (cos | sin 6 -l^-^)/2 + gin | cog ^2^-^/2 j (6 lc) 



which relate the six dynamical variables ( ) to the new variables 

(ifta, 6\, 6*2, fi, (f2)- These transformation equations are now used to evaluate the terms of 
equation fl2.10p . It is found that J2 P a * p a p — |V'o| 2 (3 + cos 6\ cos 9 2 + sin 9i sin 9 2 cos(<pi — ip 2 ))- 
Let the positions of the vortices be expressed in a coordinate-invariant manner, by defin- 
ing a unit vector to point from the centre of the vesicle towards the k th vortex. Then 
the cosine of the angle subtended by the vortices at the centre is (n x • n 2 ). In spherical 
polar coordinates, this becomes ni • n 2 = cos 6± cos 6 2 + sin 6\ sin 6 2 cos(y9i — ip 2 ). So this 
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coordinate-invariant description has arisen naturally from the model, and the first term of 
equation ( |2.10| ) becomes 

1 

a ^2 a * P a p = al^ol 2 ^ + rii ■ n 2 ). 
P =-i 



Evaluation of the other terms of equation (|2.10| ) is facilitated by noting that 

1/1 k=i 

where f is the unit vector in the direction of R. Hence 

Au 



f3pqrs a p a q a r a s 



pqrs 



57T 



ri2 2 

|4<ff! ~ h £ 

6K m=-2 



\^\ 2 Y 2 m dn 



u\ipo\ 4 
15A 



{9 (13 + 10ni • n 2 + (m • n 2 ) 2 ) - 4?r/i (3 + {n x ■ n 2 ) 2 )} . 



As one would hope, all parts of the free energy are expressible in a coordinate-invariant way. 
To complete the transformation of variables in equation Q2.9Q , a Jacobian must be evaluated. 
By calculation of a six-by-six determinant, it is established from the transformation equations 
Ul) that 

d{cL*_^, d— 1, Oq, CLq, Q^, Qi] 



l^o I (1 — ni • n 2 ) sin^i sin^ 2 . 



(6.2) 



Finally, let ipo = te la so that t is a measure of the overall real amplitude of the order 
parameter. The partition function for smectic-C vesicles becomes 

where dQk — sin 6 kd$kd(fk- As noted above, the free energy is independent of the overall 
complex phase a, which has been integrated out of the partition function, giving rise to the 
factor 2ir. Integrals over both vortex positions still remain, although it is only their relative 
separation that is relevant. Let \ be a measure of this relative separation, being equal to the 
cosine of the angle subtended by the vortices at the centre of the vesicle. Hence x ranges in 
value from -1 when the vortices are antipodal to +1 when they are coincident. This variable 
is conveniently introduced into the partition function thus: 
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Hence 



Z = 16tt 3 J d X Jltt t 5 (l - x) e~ Feff ^ (6.3) 
where F eff = A(x)t 2 + B(x)t A with A = a (3 + x) and 

B = {9(13 + 10% + X 2 ) ~~ 47r/i(3 + x 2 )} /^5uj 2 . The original set of six thermodynamic pa- 
rameters {^4, r, it, C, k} has been reduced to just three independent combinations: 

• the shifted temperature-like parameter a = r + 4imC/A, 

• the malleability /i = n 2 C 2 /Auk, 



• and the scale parameter lu = yjA/u, 

with n = 1 in this case. 

Notice that there is more phase space available to vortices at large separations than to 
vortices in close proximity; a phenomenon arising from the phase space factor (1 — x). One 
might naively predict that the thermal expectation value of X would tend to zero in the 
limit of high temperature since the vortices would spend as much time close together as 
they would in opposite hemispheres. But this is not the case. The factor (1 — X ) favours 
large inter-vortex distances, and hence (x) must tend to a finite, negative number at high 
temperature. 

The double integral of equation ( |6.3| ) is soluble; the method being given in appendix 
The solution is 

lOir^uj 2 f acu \ ( aui 

= (3 + 2 ^ a [ Wn(9-W~ W4( 27 - 2 ^)^ 

where f(x) = xe x2 eiicx and erfc is the complementary error function: erfcx = 1 — erfx. 
As x — > +oo, f(x) — > (1 — \ j1x 2 )j -Jtx and as a; — > — oo, /(x) — > 2xe x ' 2 . 

From the expressions for the partition function in equations ( |2.9| ) and ( |6.3| ), the following 
relations can be deduced: 



(6.4) 



(3 + X )t 2 ) = (E a > P ) = (J iVf^u) (6.5a) 
(13 + 10 X + X 2 )t 4 ) (6.5b) 



d\nZ 

da 
5lu 3 d\nZ 
6 duo 

from which various correlators of the variables a p are derived below. However, moments of 
X cannot be found from Z alone and, as noted in appendix [AJ, it is not easy to simply invent 
a field coupled to x to make this possible. 

Expectation values of powers of x are evaluated as follows. Returning to equation ( |6.3| ) . 
the t integral alone is easily solved, giving 



dx (6.6) 

where f'(x) is the derivative of the function f(x) defined above. Hence the moments of x 
are given by 

(0 = ^ J f_y (1 _ x)B -l/'( i ^)* £ 

which is a function of \i and (au) only. A graph of (%) against aw is plotted in figure |1| 
for the cases of the rigid sphere (/i = 0), the marginal type-I-type-II vesicle (/i = /x c ), and 
an intermediate malleability (fi = /x c /2). As increases, the vesicle becomes less rigid, and 
the lower 'kink' in the graph becomes sharper until, as fi —>■ 9/Att, (x) tends to a singular 
function for which the mean-field value (x) = — 1 is correct for negative a, as the above 
integrands are then singular for this value of x- in fact, it is generally true that mean field 
theory becomes increasingly accurate for negative aas/i approaches this critical value. The 



value /i c = 9/47T is in agreement with the critical malleability calculated in section |T|. For 
all fi, at high temperature (ie. in the large auo limit), (x) tends to 5 — 8 In 2 rs —0.545, which 
is a finite, negative number by virtue of the phase space factor discussed above. Consider 
equation (|6.6|) once more. Writing Z = J rf^exp — V(x), we see that the effective inter- vortex 
potential, 
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V(x) = - In /' I I + | In B - ln(l - X ) 



is repulsive. For this reason, (x) — > — 1 in the low-temperature limit, as this is the antipodal 
configuration. Note that V(x) is negatively dependent on /i. So the vortex- vortex repulsion 
is weakened when it is renormalized by shape fluctuations. 

Let us return to the problem of calculating expectation values from the partition func- 
tion as expressed in equation ( |6.4| ). Equations ( |6.5|) may be used to find a certain set of 
expectation values, but these are not in the form required by equations ( |4.1| ) and ( |4.2|) for 
evaluating expectation values of the vesicle's shape. How are the correlators (a* p a q ^ and 
a*a*a r a s ) to be deduced from the combinations of x an d t produced by equations ( |6.5| )? 



One can quickly convince oneself, from the form of equation (|2.10| ), with its 'momentum' 
conserving coupling, that ^OpO^ = for p ^ q. Let G P (t, x) be the mean value of a*a p , 
when thermally averaged over all configurations which have a particular value of t and X - 
This constrained average is given by 

/ a* p a p e~ FeS 5(x - Xo) S(t - t ) U q da* q da q 
p[ °' X0 ' ~ / 5( X - Xo) S(t - t ) U q da*da q ' 

Using the Jacobian of equation (|6.2|) to change the variables of integration, the constrained 
average is found to be 

G P (t ,Xo) = J a*a p 5(ni ■ n 2 - Xo)d^idQ 2 - 

And similarly, defining Q pq (t, x) to be given by thermally averaging a*a p a*a q over all states 
of a given t and X i it is found that 

G P q(to,Xo) = / Opa p a*a g 5(ni • n 2 - xo)dflidQ 2 

all other fourth-order moments being zero. Having found these constrained averages, the 
full thermal averages of a*a p and will result from averaging the functions Q p (t, x) 

and Q pq (t, x) over t and x- Hence, this is a derivation of transformation equations between 
correlators of a p and correlators of t and x- From equations ( p.5|) . 
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a\ai = t 2 (l + cos^i)(l + cos 6 2 ) 
a*_ 1 a_ 1 = t 2 (l — cos6'i)(l — cos^ 2 ) 

a^ao = t 2 (l — 2 008 6*1 cos 6*2 + rii • n 2 ). 

So, in general, the integrals required are of the form J ^(cos^i, cos 6 2 ) S(rii ■ n 2 — x)dttidQ 2 
where g{x,y) is some function. This integral then, is of a function of the positions of two 
points on a unit sphere. The points are each varied over the surface, in such a way that their 
separation remains constant. A general solution for such integrals is calculated in appendix 
O. The solution is 



g(cos 6*i, cos 9 2 ) £( n i ' n 2 ~ x)d^\dVt 2 = An ds 

Jo Jo 



VT 



where 



x — 2 s 



1 + X + V 1 ~ x) c os£ + (VI + X~ V 1 ~~ X) sin £ 



.'/ = ( V 1 + X - V 1 - X) cose + (V 1 + X + V 1 - X) sine 
Finally, the following correlation functions are found: 

(a* p a p ) = (t 2 ) + | (t 2 X 

and 



3 da 



2(al ai a* a ) = ± (13 (t 4 ) + 10 (t 4 X ) + (*Y 



u 3 d\nZ 
~9 &7~ 



uj 2 dlnZ 
2n d\x 



The derivatives are performed on equation ( |6.4| ). The order parameter is found to behave 
as follows: 



-(f\M 2 dA) = — 
u \J I au 



r 



r 



^^(27-2^) 



i(9-47T/i) 



^(27-2tt/x) 



/ / 



£(9-47T/i), 



fK 27 -2tt/^ 
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In the high-temperature limit, this tends to 3/au and, in the low-temperature limit, to 
the mean-field value of — 15aa;/(18 — 8717/). So, as p increases towards the critical value 
of p c = 9/Att, the order parameter increases in magnitude at low temperature, but not 
at high temperature, so the crossover between the two regimes becomes sharper, and 
mean-field theory becomes relatively more accurate. Figure ^ is a graph of the quantity 
(/ |^> 2 |oL4) (18 — 87r/i)/15o;, which is normalized for low temperature, against aw, for vari- 
ous values of the malleability p. 

The shape expectation values can now be found. As all of the correlators ( CLpClp\ are 
equal, substituting for (a*a g ) in equation ([O]) with (a*a p \5 pq gives 



-3 - {t\) 

\Plm) — CA / , Ipplm- 

DA; p 

But S P Jppim = 0. So all spherical harmonic amplitudes have a mean value of zero, despite 
the distorting influence of the two topological defects in ip. This is because all vortex posi- 
tions have been thermally averaged. The second moments of p are calculable, via equation 
( [4.2[ ), from the fourth moments of the order parameter field. For instance, the mean-square 
amplitude of the Y§ deformation mode is found to obey 



^(|P2o| 2 ) ^ 



7T 

+ —pau 
45 



/ 4[A( 9 - A^)]-^f{ . ™ \ - [4(27 - 2np)]-^f 



V 



15 V Wg (9-4 W Ll5V ^^(27-2^) 

r I auj \ r I olj \ 

Wb( 9 -m/ \ 75(27=2^; 

The constant ^ is due to thermal excitation, and the other term arises from deformation by 
the in-plane order. As p — > p c , this function goes to zero for positive a, and to 12 g^ (jf^Jl) 
for negative a. The latter expression tends to infinity as expected, since deformations become 
large as the marginal type-I-type-II case is approached. The function (p c — p) 2 n {\p2o\ 2 ), 
which is normalized for low temperature, is plotted against au in figure [| for the cases 
p = p c /2 and p = p c . 

Notice that all of the measurable quantities calculated in this section are functions of 
p and the combination aw. Hence, temperature, measured from the mean-field transition, 
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scales asw 1 . So 'transitions' from low- to high-temperature behaviour become sharper for 
larger vesicles (ie. when A is large compared with u). 



VII. HARTREE-FOCK APPROXIMATION 

Up to this point, the only approximation employed, other than the small amplitude ap- 
proximations implicit in the Ginzburg-Landau model, has been the confinement of the order 
parameter field ip to a 2(2n + l)-dimensional phase space of the lowest eigenfunctions of 
the gradient operator. A further approximation is now introduced, in order to find expecta- 
tion values of the n-atic order-parameter field for all of the values of n under consideration. 
The Hartree-Fock method will be used to produce an infinite, but incomplete, perturbation 
expansion, valid at high temperature - the opposite limit to that studied in |TT] . 



Firstly, the bare, or Gaussian propagator h is calculated from the quadratic Hamiltonian 

H = a* p a P 
p 

and is found to be independent of p. 

v / a*a q e- Ho \[ r da* r da r 5 pq , 

aa n ) = r „, — : : = = nn h 



P^/o~ Je-^u r da* r da r " a ~ m p 
Writing the free energy as F = H. +li,i where Hi = J2 pq rs P" s 

the partition function 

becomes 

where Z is the partition function for the quadratic Hamiltonian and () indicates a thermal 
average with respect to this Hamiltonian. Hence, the renormalized propagator (ie. the full 
thermal average, can be written as 



hp = (a* p a p ^ 



ct r , el* n a 



Taylor-expanding the exponentials and applying Wick's theorem, the terms in this formula 
may be represented by connected, two-leg graphs of all possible topologies, in which di- 
rected lines represent the bare propagators with 'momentum' p, and the four-point vertices 

23 



each carry a factor of /3 eff . Momentum is conserved at the vertices, since f3p B qrs vanishes if 
p + q 7^ r + s. The numerator is the sum of all connected and disconnected two-leg graphs, 
but the disconnected parts cancel with the denominator, which is the sum of all zero-leg 
graphs. This explanation is brief, since the Feynman- diagram expansion is given in many 
standard texts ]TS| . 

Truncating the expansion at one-loop graphs yields 

connected 

ht ~ h t — ^2 Ppqrs (at a p a q a r a sQ>t} 
pqrs 

connected 

= h t — ^2 Ppqrshth p h g (5 tr + 8 ts )(8 pr 5 qs + 5 ps 5 qr ) 

pqrs 

represented by 




> fa — =» + 

t t 



where a double line represents a renormalized propagator. 

The Hartree-Fock scheme renormalizes the bare propagator with not just one loop dia- 
gram, but an infinity of diagrams belonging to a certain topological family. In particular, 
all those diagrams containing loops which do not span a vertex. This approximation to the 
renormalized propagator is defined implicitly by the equation 




^> ^ _|_ 

Y~ , 

p jTv 

which can be thought of as iteratively replacing every bare propagator with the one-loop 



correction. The scheme generates graphs with the correct counting factors IT]]. This 
pictorial equation corresponds to 
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hp ~h p (l- 2h p J2Wf qqP + PZjOj 



a set of (2n + 1) second order coupled polynomial equations. These generate many spurious 
solutions. But it is observed, from calculation of (3 p ° qrs that J2 q {P pqqp + Pp%q) is independent 
of p, and hence there is always a solution for which the (2n + 1) quantities h p are equal, 
and positive. Henceforth, let h p = h Vp. This physical solution is given by the quadratic 
equation 



h 2 ( 2 E0%p + + ah -1=0 Vp. 



(7.1) 



The coefficient of h 2 is now calculated. 
Firstly, 



E /%, = ^ M E W^-^-Jr^ E E 1 j /<WTd« /*r E 

The function YZ= _ n I0g| 2 5 which appears twice in the above equation, may be written as 



(2n + l)! ^ sin 2 - (f) cos 2 ( 2 — ) (f) 
4vr ^ r!(2n-r)! 

which is in the form of a binomial expansion of (sin 2 (|) + cos 2 (|)). It is therefore just a con- 
stant; J2 q |0g| 2 — (2w + l)/47r; and hence orthogonal to all non-s-wave spherical harmonics. 



It follows that 



E fipqqp ( n + 2 ' 



A 

Secondly, the sum J2q Ppqpq * s found. This also has two parts, deriving from equation 
Q2.12| ), and the first has already been calculated, since (3 pqpq = (3 pqq p. However, the order of 
indices is relevant to the second part, which does not vanish this time. There is no obvious 
method for simplifying this part, so the integrals have been calculated as required for each 
value of n of interest, with the result that 

E« 9P + «,) = (2n + 1)^ - 8j ^f(n) 
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where /(l) = i, /(2) = ^, /(4) = jgg and /(6) = ^gfg. So the solution to equation 
(0) ^ 



2co e V £ 2 
where 



1 Vp. (7.2) 



e = 4Jn + ~ — 167r/(n)//. 



As before, /x is the 'malleability', n 2 C 2 / Auk, and is the 'scale parameter', ^jA/u. From 
the definitions, the magnitude of the order parameter is 



/" \^j\ 2 dA = Y J K = (2n+l)h. 



In figure £|, a graph oieh/uj against au/e is plotted in bold. Recall that a is a temperature- 
like variable, measured with respect to the shifted mean-field transition temperature r c . 
The figure shows that, within the lowest Landau level, and Hartree-Fock approximation 
schemes, the transition is removed. Instead there is a gradual change from a low-temperature 
region with a high degree of order, to a poorly-ordered high temperature region, and this 
picture may be closer to the truth. Of course, the Hartree-Fock approximation should only 
be trusted at high temperature (au/e ^> 1), as it is rigorously correct to 0{u). At low 
temperature, equation (|7.2|) tends asymptotically to h — —aoj 2 /{n + | — 167r/(n)/x), which 
is at variance with the correct mean-field expression h — > — a/6/?oooo- -^ or n = 1> ^ ne f° rmer 
expression becomes — 3aa; 2 /(18 — 167171) and the latter — 5au 2 /(18 — 8717/). Nonetheless, 
the Hartree-Fock expression is a useful qualitative indicator of the system's behaviour. For 
comparison, the correct, lowest-Landau- level solution for n — 1, as calculated in section 0, 
is plotted as dashed lines in figure |], for various values of the malleability fi. Note that 
some /i-dependance is contained in the function e, which scales the axes of the figure. It 
is apparent from the graphs that the Hartree-Fock solution is correct at high temperature, 
and fairly poor at low temperature, except for the special case when /1 = 9/147T, for which 
the asymptotic gradient is in agreement with the mean-field value. 
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The ordinate of figure (| is a measure of the order parameter, under some normalization, 
while the abscissae is the temperature parameter a in units of e/u. Hence it is immediately 
apparent for the expression for u that the crossover between high- and low-temperature 
regimes is sharper on larger vesicles. The crossover apparently becomes infinitely sharp 
at a finite value of the malleability, when /i = /i c = (2n + l)/(327r/(n)). Once again, it is 
apparent, from reference to the correct values of \x c calculated in section fH, that the Hartree- 
Fock approximation gives qualitatively correct, but quantitatively inaccurate results. 

The shape expectation values can now be found. Substituting for (a*a q ^ in equation 
Q4.1D with hSpq gives 

(Pirn) — r, a Ipplm- 
ZiX l p 

But J2plppim = 0. So all spherical harmonic amplitudes have a mean value of zero, de- 
spite the distorting influence of the 2n topological defects in %p. This is because all vor- 
tex positions have been thermally averaged. The second moments are found, via equa- 
tion ( |4.2[ ), from the fourth moments of the order parameter field. From Wick's theorem, 
a*a*a r a s ^ = (5 pr 5 qs + 5 ps 5 qr )h 2 . So equation ( |4.2j ) becomes 



Pl imi Pl2 

For instance, the auto-correlation function (IP20I 2 ) for n = 1 is evaluated as 

|2\ _ 1 ( A1X ^12 



Hence deformations increase with the magnitude of the order parameter, as expected. 



VIII. CONCLUSION 

A new, simple expression has been derived for the coupling between the lowest Lan- 
dau levels of intrinsic n-atic order and the shape of a genus-zero membrane. It has been 
shown, within the lowest-Landau-level approximation scheme, that all spherical harmonic 
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amplitudes of deformations of such a vesicle have zero thermal average, but that their mean 
squares have a contribution from deformation by the order, and a contribution from thermal 
excitation. Only a finite number of spherical harmonics are coupled to the lowest Landau 
levels of order. It has been established that this order is capable of fundamentally alter- 
ing, and perhaps bursting a vesicle of finite rigidity. This occurs at a critical value of the 
'malleability', for which the model is strongly analogous to that of a marginal type-I-type-II 
superconductor. As this critical value is approached, mean-field theory becomes increasingly 
accurate. The model for fluid vesicles with vector order has been solved exactly, using no 
further approximations, and the general result is established that the crossover from high- 
to low- temperature behaviour becomes sharper for larger vesicles. The solution also con- 
firms that mean-field theory becomes increasingly accurate as the critical malleability is 
approached. The counter-intuitive result is ascertained that the two topological defects in 
the order are more likely to be found far apart than close together, even at high tempera- 
ture, as the phase space volume element increases with defect separation. The Hartree-Fock 
approximation is valid at high temperature and significantly inaccurate at low temperature, 
but exhibits qualitatively correct behaviour. 
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APPENDIX A: SOLUTION OF THE PARTITION FUNCTION DOUBLE 

INTEGRAL 

The integral required in equation (|6.3|) is 
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where A = a(x + 3) and B = b(x 2 + 10x + 13) — c(% 2 + 3). Expanding the first exponential 
as a sum and using t 2m+5 e~ Bt4 dt = T( m ^-)/AB m i 1 gives 



m=0 m - ^ 



i3 2 



Using the indefinite integral 



i 



B n+1 



-(3 + X ) 2n . , 

rfv = — — — — h const. 

A 2n(26 + 3c)B n 



gives 



IE 



-a) m (^)! 



'„ (m+ l)!(26 + 3c) 



(6_ c )^ V66-c 



m + l 



The exact form of the function -B(x) is crucial to the solubility of the indefinite integral above. 
In particular, it is not easy to invent new fields coupled to the dynamical variables to aid com- 
putation of their expectation values. Now, applying the identity (y)!/m! = ^/2 m ( !n ^)\ 
allows the infinite sum to be separated into the standard forms of exponentials and confluent 
hypergeometric functions thus: 



1 



16(26 + 3c) 
+ 



2^ 



la 



exp 



4(6 - c) 



1F1 1, 



4y^F 

2a 



exp 



6b — c 



1F1 1, 



Sb-cJ X ~ X \~'2 J 6b- c) \b-cj \ 2 4(6 -c) / 
The identities i-F\(a, c; x) = e^iF^c — a, c; — x) and erf (x) = ^7§i-F±(§, |; — x 2 ) are used |T9 



to give the result 



7T 



4a(26 + 3c) 



x e x erfc x 



where erfc is the complementary error function: erfcx = 1 — erfx. 



APPENDIX B: SOLUTION OF THE CONSTRAINED, TWO-POINT, 
SPHERICAL SURFACE INTEGRAL 



The integrals required for thermally averaging combinations of a v over states of a given 
t and x are °f the form 
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J = J g(cos 0i, cos 6 2 ) S(ni ■ n 2 — x)d^\dVl 2 

where g(x,y) is some function. The integral is over the positions of two points on a unit 
sphere. The unit vectors n\ and n 2 point from the centre of the sphere to the points, each 
of which is varied over the surface, in such a way that rt\ • n 2 has the constant value \- I* 1 
spherical polar coordinates, the solid angle elements are dVL k = d(cos9k)dipk- Hence 

J = J d(cos#i) J d(cos9 2 ) <?(cos#i, cos# 2 ) 

/ dip\ / d(<fi — <f 2) 5(cos^i cos ^2 + sin 6*i sin# 2 cos(y?i — if 2) — x) 
Jo Jo 

and 

/■27T 

/ dip 5 (cos 61 cos 9 2 + sin 9\ sin # 2 c °s <f — x) 
Jo 

n f 1 d(cosip) „ . ./,./, \ 

= 2 / (cos 1 cos # 2 + sin 9 ± sin # 2 cos </? — x ) 

7-i sin 6*i sin 6*2 sine/? 

26 f ( x r cose l C o S e 2 y 

y \ sin Oi sin #2 / 



\j sin 2 #1 sin 2 9 2 — (x ~ cos #1 cos $2 
where 0(x) is the Heaviside step function; 0(a; < 0) = 0, Q(x > 0) = 1. Now, setting 
x = cos Q\ and y = cos 9 2 gives 

J = 4n[ 1 dx[ 1 dy , 9 =M Q< {X ~ ^ ^ 



1 J-i 



^(l-x 2 )(l-y 2 )-( X -xy) 2 \(l ~ x 2 )(l - y 2 ) ) ' 

The region of the x-y plane for which the Heaviside function has a positive argument is the 
interior of an ellipse, contained within the unit square of integration. The transformation 



2x = (y/l + X + V 1 " X)u + (y/l + x - y/l - X> 
2y = (y/TT^ - yjl-x)u + (yfi+X + ^ - X> 



with Jacobian 



d(x,y) 



X 2 



d(u, v) 

maps the ellipse onto a unit circle in the u-v plane thus: 
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-VT^ \J\-u 2 - 



Hence a final transformation to polar coordinates (s, £) gives 



J = 4tt 



ds dC 
Jo Jo 



■g(x,y) 



with 



x — 2 s 



y=r s 



1 + X + V 1 ~ X) cos £ + (y 1 + X ~ V 1 - X) si n C 



So, for instance, 



Js(n 1 ■ n 2 — x) dQ\dQ:2 = 8n 2 



and 



J 3os 9i cos 2 5(rii • n 2 — x) dQidQ 2 = 



8tt 2 X 
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FIGURES 

FIG. 1. Graph of the expectation value of the cosine of the angle subtended at the vesicle's 
centre by the two vortices (x) against the temperature-like combination qw, for the cases \x = 
(the rigid sphere), /i = /i c (critical malleability), and fj, = fi c /2. The vortices are antipodal at low 
temperature, and tend to be well separated even at high temperature due to the shape of the phase 
space. 

FIG. 2. Graph of the normalized order parameter (J \tp 2 \dA) against the temper- 

ature-like quantity aw, for various values of the malleability [i. As \x approaches the critical 
value /i c = 9/47T, low-temperature order increases, since the vesicle is intrinsically flattened, but 
high-temperature order is independent of \i. Hence mean-field theory becomes accurate in the 
marginal type-II-type-I limit, [i = [i c . 

FIG. 3. Graph of the mean-square amplitude of the mode of deformation multiplied by 
{n — Hc) 2 k, against the temperature-like quantity auj, for the cases fi = /x c /2 and fi = \i c . 

FIG. 4. Graph of eh/u against au/e, showing the non-singular progression from the 
low-temperature, highly ordered state to the high temperature poorly ordered state. The 
Hartree-Fock solution is show as a continuous, bold line. For comparison, the more accurate 



solution for n = 1, as calculated non-perturbatively in section VI, is shown dashed, for various 
values of the malleability [i. 
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